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ABSTRACT 

We infer the UV luminosities of Local Group galaxies at early cosmic times (z ~ 2 
and z ~ 7) by combining stellar population synthesis modeling with star formation 
histories derived from deep color-magnitude diagrams constructed from Huhble Space 
Telescope {HST) observations. Our analysis provides a basis for understanding high-z 
galaxies - including those that may be unobservable even with the James Wehh Space 
Telescope {JWST) - in the context of familiar, well-studied objects in the very low-z 
Universe. We find that, at the epoch of reionization, all Local Group dwarfs were 
less luminous than the faintest galaxies detectable in deep HST observations of blank 
fields. We predict that JWST will observe z ~ 7 progenitors of galaxies similar to 
the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may 
already be observing such galaxies, highlighting the power of gravitational leasing. 
Consensus reionization models require an extrapolation of the observed blank-field 
luminosity function at z « 7 by at least two orders of magnitude in order to maintain 
reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf 
spheroidal galaxies to be contributors to the ionizing background at z ~ 7. Combined 
with numerical simulations, our results argue for a break in the UV luminosity function 
from a faint-end slope of a ^ —2 at Muv ^ —13 to a ^ —1.2 at lower luminosities. 
Applied to photometric samples at lower redshifts, our analysis suggests that HST 
observations in lensing fields at z ^ 2 are capable of probing galaxies with luminosities 
comparable to the expected progenitor of Fornax. 

Key words: Local Group - galaxies: evolution - galaxies: high-redshift - cosmology: 
theory 


1 INTRODUCTION 

Faint galaxies below the detection limits of current observa¬ 
tories are necessary contributors to the ionizing background 
in the galaxy-dominated models of cosmic reionization that 
are cnrrently favored. Deep-held observations with the Hub¬ 
ble Space Telescope {HST) only reach « 0.11/*(z ~ 7), 
while essentially all models of reionization require galax¬ 
ies 10 — 1000 times fainter to contribute to the ionizing 
background (Alvarez et al. 2012; Kuhlen & Faucher-Giguere 
2012; Duffy et al. 2014; Robertson et al. 2015; Bouwens et al. 
2015a). Probing these likely drivers of reionization is a prime 
motivation of the James Webb Space Telescope {JWST), yet 
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even JWST is unlikely to reach the faintest (and most nu¬ 
merous) of these sources (see Section 3.2). 

Resolved-star observations of Local Group galaxies pro¬ 
vide an alternate path for learning about the faintest galax¬ 
ies at the epoch of reionization, via “archeological” studies 
of their descendants. Local galaxies also provide important 
benchmarks against which data at a variety of redshifts can 
be compared. Accordingly, observations in the Local Group 
have long informed our understanding of faint galaxies at 
early times (Bullock et al. 2000; Freeman & Bland-Hawthorn 
2002; Ricotti & Gnedin 2005; Madau et al. 2008; Bovill & Ri- 
cotti 2011; Brown et al. 2012; Benftez-Llambay et al. 2015). 

The relation between high-z star formation rates 
(SFRs) and dark matter halo masses for the progenitors 
of Local Group dwarfs is of particular interest for under¬ 
standing the reionization epoch. Boylan-Kolchin et al. (2014, 
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hereafter B14) showed that the requisite inefficient star for¬ 
mation in low-mass dark matter halos at 2 = 0 proves diffi¬ 
cult to reconcile with expected physics of the high-redshift 
Universe in general and reionization models in particular, 
as these often require relatively efficient star formation in 
low-mass halos at early times. Indeed, Madau et al. (2014) 
suggested that the reionization-era progenitors of present- 
day dwarf galaxies may have been among the most efficient 
sites of conversion of gas into stars on galactic scales. Effi¬ 
cient galaxy formation at high redshift has low-redshift im¬ 
plications, however: based on counts of remnants of atomic 
cooling halos from 2 = 8 that survive to 2 = 0 in simulated 
Local Groups, B14 argued that the UV luminosity function 
(LF) should break to a shallower value than the observed 
a ~ —2 for galaxies fainter than Muv ~ —14 (correspond¬ 
ing to halo masses of Mvir 10® Mq). 

Weisz et al. (2014b, hereafter W14) further demon¬ 
strated the power of combining local observations with stel¬ 
lar population synthesis models in confronting LFs in the 
high-redshift Universe: they showed that one can construct 
LFs that reach several orders of magnitude fainter than is 
currently observable. The results of W14 indicate that UV 
LFs continue without a sharp truncation to much fainter 
galaxies [M\jy{z ~ 5) « —5], showing how near-held obser¬ 
vations can inform our understanding of galaxies that are 
not directly observable at present. 

In this paper, we combine the approaches of B14 and 
W14. Using resolved star formation histories (SFHs), we 
model the UV luminosities that Local Group galaxies had 
at earlier epochs. We compare these to existing HST data 
at 2 ~ 2 (from UV dropout galaxies in blank helds and 
cluster lensing helds) and at 2 ~ 7 (from the Ultra-Deep 
Field (HUDF) and the Frontier Fields). As we show below, 
this approach holds the promise of both placing the faintest 
observable galaxies at high redshift into a more familiar con¬ 
text and understanding the possible role of well-studied Lo¬ 
cal Group galaxies in high-redshift processes such as cosmic 
reionization. 

2 MODELING LOCAL GROUP GALAXIES AT 

HIGHER REDSHIFTS 

2.1 Deriving Muv(z) 

To interpret our observations of Local Group galaxies in a 
high-redshift context, we follow the methodology described 
in W14. This procedure is summarized here; for further de¬ 
tails, see W14 and references therein. 

We begin with SFHs derived from resolved star color- 
magnitude diagrams (CMDs) constructed from HST imag¬ 
ing. The majority of the SFHs we use are based on analysis 
of archival HST /WFPC2 imaging; these data, the SFHs, 
and details on how the SFHs were computed are presented 
in Weisz et al. (2013) and Weisz et al. (2014a). We select 
Local Group galaxies that span a wide range in present-day 
stellar mass (2 x 10® — 2 x 10® Mq) and have HST-hased 
CMDs that extend below the oldest main sequence turn-off 
(MSTO), which enables a precise constraint on the stellar 
mass at all epochs back to 2 = 5. The archival WFPC2 data 
meet the MSTO criteria for galaxies within ~ 400 kpc of 
the MW, and of these, we select a representative subset. We 


supplement this dataset with SFHs of Leo A (presented in 
Cole et al. 2007) and IC 1613 (presented in Skillman et al. 
2014), which are based on newer imaging that includes the 
ancient MSTO. Our sample includes SFHs for all known 
bright {Lv > 10® Lv,©) MW satellite galaxies with the ex¬ 
ception of Sextans (for which there is no sufficiently deep 
HST imaging), as well as several Local Group dwarf irreg¬ 
ular galaxies. SFHs from Local Group galaxies are “non- 
parametric” in the sense that they are the sum of simple 
stellar populations (e.g., Dolphin 2002) as opposed to im¬ 
posed analytic prescriptions (e.g., r models). 

The SFH of each galaxy is input into the Flexible Stel¬ 
lar Population Synthesis code (Conroy et al. 2009; Conroy 
& Gunn 2010) in order to generate UV and V-band flux 
profiles over time. In both of these steps, we use the Padova 
stellar evolution models (Girardi et al. 2002, 2010) and a 
Kroupa (2001) initial mass function (IMF). We assume a 
constant metallicity of 0.2 Z© and no internal dust extinc¬ 
tion (see the discussion in Section 2.3). The latter assump¬ 
tion appears to be in broad agreement with the apparent 
dust-free spectral energy distributions of faint, high-redshift 
galaxies (e.g., Bouwens et al. 2012; Dunlop et al. 2013) and 
with expectations from simulations (Salvaterra et al. 2011). 
Finally, in order to generate the simulated fluxes as a func¬ 
tion of redshift, we require the simulated V-band luminosity 
at 2 = 0 to match observations of Local Group galaxies as 
listed in McConnachie (2012). This normalization accounts 
for the fact that the HST field of view does not always cover 
the entire spatial extent of a nearby galaxy. 

We initially generate predicted fluxes assuming the fidu¬ 
cial SFH, i.e., a constant SFH over each time bin. However, 
as suggested by a variety of simulations® (Stinson et al. 2007; 
Ricotti et al. 2008; Governato et al. 2012; Zolotov et al. 2012; 
Teyssier et al. 2013; Dominguez et al. 2014; Power et al. 
2014; Oiiorbe et al. 2015) and observations (e.g., van der 
Wei et al. 2011; Kauffmann 2014), the SFHs of low-mass 
galaxies at high redshift are likely to fluctuate on timescales 
(~ 10-100 Myr) that are shorter than what is directly avail¬ 
able from the fossil record (which provides a time resolution 
of ~ 10-15% of a given lookback time for CMDs that include 
the oldest MSTO). 

To account for short duration bursts of star formation, 
we have modified the fiducial SFHs of our Local Group sam¬ 
ple to include short duration (~ 10-100 Myr) bursts. That is, 
in each time bin, we insert bursts with a specified amplitude 
and duration that are stochastically spaced in time, with 
the requirement that the total stellar mass formed in each 
time bin matches that of the fiducial SFH. For the purpose 
of this exercise, we have adopted two representative burst 
models: 200 Myr duration with amplification factor of 5 and 
20 Myr duration with amplification factor of 20. In both 
cases, we assume that 80% of the star formation occurs in 
the bursting phase. Each modified SFH is therefore a series 
of 20 (200) Myr bursts with an amplitude of 20 (5) times 
the average SFH instead of the fiducial constant SFH in a 


^ Not all simulations agree on this point: for example, Vogels- 
berger et al. (2014) were able to match many observable proper¬ 
ties of dwarf galaxies with galaxy formation models that result in 
much smoother star formation rates as a function of time. 
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redshift redshift redshift 

20 10 7 6 5 20 10 7 6 5 20 10 7 6 5 



13.6 13.4 13.2 13.0 12.8 12.6 13.6 13.4 13.2 13.0 12.8 12.6 13.6 13.4 13.2 13.0 12.8 12.6 

Lookback Time [Gyr] Lookback Time [Gyr] Lookback Time [Gyr] 


Figure 1. Star formation rates calculated in bins of 200 Myr (light gray histograms) and 20 Myr (dark gray histograms) for simulated 
galaxies from Onorbe et al. (2015; left) and Fitts et al. (in preparation; center and right) up to 2 : = 5. At the present day, each galaxy 
has 10® < M*(z = 0)/Mq < 5 x 10® and is hosted by a halo with = 0) = lO'^® Af©- The dotted horizontal line shows the mean 

SFR over the period plotted. Averaged over 200 Myr periods, the SFRs appear to be increasing to 2 = 7. On 20 Myr timescales, they 
are much burstier and fluctuate strongly. These simulated SFHs motivate the burst parametrizations we use in our modeling. 


given time bin (as noted above, each time bin has a width 
of ~ 10-15% of the corresponding lookback time). 

The choice of these parameters is guided by cosmolog¬ 
ical simulations. Figure 1 shows examples of such episodic 
SFHs for three simulated dwarf galaxies at 2 > 5. Each 
was run using Gizmo (Hopkins 2015) with meshless finite- 
mass hydrodynamics at very high force and mass resolution 
and using the FIRE implementations of galaxy formation and 
stellar feedback (Hopkins et al. 2014). The dwarf in the left 
panel is a version of the “dwarf-early” galaxy presented in 
Onorbe et al. (2015) and Wheeler et al. (2015), while the 
dwarfs in the center and right panels are of halos with nearly 
identical masses selected from a 35 Mpc volume and resimu¬ 
lated at very high resolution (from A. Fitts et ah, in prepa¬ 
ration). The average SFRs of these galaxies at 2 > 5 are 
(0.1 — 1) X 10“® Mq yr“^ and their stellar masses at 2 = 0 
are (1.8-5) x 10® M©. 

Figure 1 shows the early (2 > 5) star formation 
in each halo, with the light gray (dark gray) histogram 
showing the star formation averaged over 200 (20) Myr 
periods. While the variation when averaged over 200 Myr 
periods is typically a factor of 2-5, the variation on 20 
Myr timescales is frequently a factor of 10-20 and can 
even exceed 100. This is further corroborated by SFRs for 
a variety of halos over the range 11 > 2 > 5 presented 
in Ma et al. (2015). While the masses of their halos 
[Mvir (2 = 6) ~ 10^®“^®'® M©] are generally larger than 
those considered here [Mvir (2 = 6) ~ 10®'®“®'®], the SFHs 
are qualitatively similar in their very bursty, episodic nature. 


2.2 From Muv(z) to P(Muv) 

The result of this modeling is a distribution of magnitudes as 
a function of lookback time, Muv(t), or redshift, M\jy{z), 
for each galaxy. The probability distribution P(Muv) for 


each galaxy is then given by P(Muv) oc Muv(«) P{z), i.e., 
it is a distribution over the modeled UV flux as a function of 
time, weighted by the probability that the galaxy falls into 
a sample with redshift selection function P{z). The weights 
can correspond to the photometric redshift distribution from 
an observational sample, to some modihcation thereof, or to 
an arbitrary selection function. 

In what follows, we adopt a Gaussian centered at 2 = 7 
with width = 0.475 as the redshift selection function 
P(z) for our 2 ~ 7 galaxies. This produces a good match to 
the photometric redshift distribution from Finkelstein et al. 
(2014) (after excluding the secondary peak at 2 ~ 1.4 origi¬ 
nating from detections of a 4000 A break rather than the 
Lyman break). Similarly, we use a Gaussian centered at 
2 = 1.91 with a width of CTz =0.21 as the redshift selection 
function P{z) for our 2 ~ 2 sample, as this reproduces the 
photometric redshift distribution from Oesch et al. (2010). 
At both 2 ~ 7 and 2 ~ 2, the width of P{z) is smaller than 
the time resolution of the observed SFHs. 

Figure 2 shows the probability distribution P(Muv) for 
the Fornax dwarf spheroidal (dSph) galaxy, computed from 
1000 realizations of its SFH, at 2 ~ 7. The left panel shows 
the case of 20 Myr bursts while the right panel shows the 
result for 200 Myr bursts. The distributions are shown as 
shaded histograms; the shading gives the value of the cumu¬ 
lative distribution function P(< Muv) for each Muv- Taking 
the instantaneous value of Muv at 2 = 7 in each realization 
results in a cumulative distribution that is very similar to 
that shown here in each case. 

While the full range of Muv spanned by the distribu¬ 
tions is set by the burst fluctuation amplitude and is there¬ 
fore roughly the same in the case of 20 or 200 Myr bursts, the 
distributions themselves differ substantially. 20 Myr bursts 
result in a bimodal distribution with a prominent peak ~ 0.5 
magnitudes from the faint end of the distribution and a 
secondary peak ~ 1.0 magnitudes from the bright end of 
the distribution. The median value falls at ~ 1 magnitude 
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Figure 2. Probability distribution for Myy of Fornax, given the photometric redshift distribution of Finkelstein et al. and assuming that 
Fornax’s high-redshift star formation occurred in bursts with a characteristic timescale of 20 Myr (left) or 200 Myr (right). The colors 
indicate the cumulative probability distribution in each panel. In the case of 20 Myr bursts, the frequent and relatively high-amplitude 
bursts result in a bimodal probability distribution function (with peaks corresponding to the burst and inter-burst periods). The 200 
Myr bursts, which are modeled with smaller burst amplitudes, result in a distribution that is less broad and is unimodal. This is a 
direct consequence of our smaller assumed burst amplitude for 200 Myr bursts; the resulting SFH is therefore much closer to the fiducial 
constant SFH. Taking the instantaneous value of Muv at z = 7 in each realization results in a cumulative distribution that is very 
similar to that shown here in each case. 


from the faint end of the distribution. In the case of 200 
Myr bursts, the distribution is noticeably different: there is 
a much stronger peak near the median of the distribution, 
with relatively high probabilities of being up to half a magni¬ 
tude brighter or a magnitude fainter than the median value 
of Muv- The net result is that the median probability in 
the case of 200 Myr bursts is approximately 0.5 magnitudes 
brighter than in the (likely more realistic) case of 20 Myr 
bursts (e.g., Weisz et al. 2012). This is similar to the value 
of Muv that would be obtained by assuming a constant SFR 
over the entire period to z = 5, i.e., the parameters we adopt 
for 200 Myr bursts act mainly to distribute Muv about the 
median value for the constant SFH. 

Given the expectations from observations and cosmo¬ 
logical simulations, we adopt 20 Myr bursts as our default 
model for star formation and compute UV magnitudes from 
the resulting probability distribution derived from 100 inde¬ 
pendent realizations per galaxy. The results at z ~ 2 and 
z ~ 7 are listed in Table 1. We quote median values and 
confidence intervals comprising the minimum range in Muv 
containing 68% and 95% of the cumulative probability. Since 
virtually all of the galaxies have P(Muv) distributions that 
look like Figure 2, these ranges are usually asymmetric about 
the median. 

An immediate connection that we can make between 
high redshifts and current-day quantities is to compare ob¬ 
served V-hand magnitudes at z = 0 with modeled values 
of M-[jy{z = 7). Such a comparison is shown in Figure 3. 
Perhaps unsurprisingly, there is a clear correlation between 
high-redshift UV luminosity and present-day U-band lumi¬ 
nosity for the dwarf galaxies studied here. The slope of this 
correlation is 0.71, indicating that the high-redshift UV lu¬ 
minosities span a narrower range than the z = 0 U-band 


luminosities. While there is non-negligible scatter in this re¬ 
lation, which is related to the diverse early SFHs of the 
low-mass dwarfs in particular (e.g., figure 12 of Weisz et al. 
2014a), the correlation is well-defined over the range of lu¬ 
minosities studied here. We caution against extrapolating 
this relation to significantly higher luminosities, as it is un¬ 
likely to remain linear indefinitely. The balance between 
M\jy{z ~ 7) and Mv{z = 0) is set by the stellar mass 
formed at early times versus over a galaxy’s lifetime, and 
very massive galaxies today have predominantly ancient stel¬ 
lar populations (Thomas et al. 2005), which is more similar 
to low-mass classical dSphs in the MW than to galaxies at 
intermediate mass (such as the MW itself). 

For our sample of galaxies as a whole, Muv(z = 2) 
is reasonably well-correlated with Muy{z = 7) and with 
Mv{z = 0). There are notable exceptions, however, includ¬ 
ing Leo A, Leo I, Leo T, and Sagittarius. These galaxies show 
much lower UV luminosities at z ~ 2 than would be naively 
expected based on their present-day luminosities. The origin 
of this difference lies in the SFHs of these galaxies: while all 
of the galaxies in our sample have ancient (z ~ 7) star forma¬ 
tion, the galaxies listed above appear to have gone through 
relatively quiescent phases at lower redshifts (see figure 7 
of Weisz et al. 2014a). These periods of quiescence result in 
dramatically reduced UV luminosities at those times. 

As discussed at the beginning of this section, we have 
only included a representative set of Local Group galaxies. 
They are mainly selected based on diversity in stellar mass 
and morphology at z = 0, and all have the requirement of 
reaching below the oldest MSTO in order to provide the best 
possible constraints on the stellar mass formed by z = 5. Fu¬ 
ture work analyzing the GMDs of newer, deep HST obser- 
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Table 1. Observed V-band magnitudes of Local Group dwarfs 
(column 2), along with modeled UV magnitudes at 2 ~ 2 (column 
3) and 2 ~ 7 (column 4). The Muv values quoted in columns 
3 and 4 are the medians of the probability distributions, while 
the errors are the minimum range containing 68 % (95%) of the 
cumulative probability for ALpy- 


Name 

My (2 = 0 ) 

Mjjy{z 2) 

Muv (2 ~ 7) 

Carina 

-9.1 ± 0.5 

„+0.6{2.4) 

( 2 . 0 ) 

. q+0.9(1.3) 
°-^-0.9(1.5) 

CVn I 

- 8.6 ± 0.2 

^ 0+3.9 (3.9) 
^•^-1.4(3. 2 ) 

q .+0.9 (0.9) 

^•'^-o.e ( 1 . 6 ) 

Draco 

- 8.8 ±0.3 

q .+0.6 (2.3) 

^ -1.7 (2.2) 

Q 9+0.9 (0.9) 
^•^-0.5 (1.4) 

Fornax 

-13.4 ±0.3 

-,9 ,+0.9(0.9) 
-^^•-^- 0.6 ( 1 . 6 ) 

^^•'^-0.5 (1.4) 

IC 1613 

-15.2 ± 0.2 

,9 .+ 0 . 6 ( 1 . 6 ) 

-^^•’^- 1.2 ( 1 . 8 ) 

,+0.9(1.0) 

J-'^-^_0.7(1.6) 

Leo A 

- 12.1 ± 0.2 

. 7+0.7 (0.9) 
'^•'-0.5(1.2) 

,q .+ 1 . 0 ( 1 . 0 ) 

^^•'^-0.6 (1.3) 

Leo I 

-12.0 ±0.3 

. ,+0.9 (0.9) 
'^•-^-2.2 (4.5) 

,+l. 0 ( 1 . 0 ) 
^-^•-^-0.6 (1.4) 

Leo II 

-9.8 ±0.3 

^ 7+1.2 (1.5) 
®’'- 0 . 8 ( 2 . 2 ) 

„ „-l-0.7(0.7) 
“■®-0.5 (1.5) 

Leo T 

-8.0 ±0.5 

, .+0.4 (0.8) 
^•^-0.4(0. 8 ) 

7q+0.9(0.9) 

‘ ‘^-0.6 (1.5) 

LMC 

-18.1 ± 0.1 

1 '5 {0-9) 

-^'^•^- 0.6 ( 1 . 6 ) 

_,.g+ 1 . 0 (l.l) 

^^■®-0.6(1.4) 

Phoenix 

-9.9 ±0.4 

qq+ 1 . 0 ( 1 . 0 ) 

^•'^-0.6 (1.5) 

1 n 

^^•^-0.6 (1.5) 

Sagittarius 

-13.5 ± 0.3 

+0.4 (0.8) 
^•^-0.4 ( 0 . 8 ) 

,9 o+0.9(0.9) 
^^•'^-0.5 (1.5) 

Sculptor 

-11.1 ± 0.5 

+0.5 (2.3) 

-^^•^-1.5 (1.9) 

q + 0.6(1.3) 
-*--*■•^-1.2 (1.3) 

SMC 

-16.8 ± 0.2 

+ 0 . 8 ( 1 .3) 

-^^•'^- 1.0 ( 1 . 8 ) 

_ +0.4(1.3) 

-‘-^•-‘■- 1.2 ( 1 . 2 ) 

Ursa Minor 

- 8.8 ±0.5 

+ 1.7(1. 6 ) 

'^•^-2.9 (5.1) 

„ q+ 1 . 0 ( 1 . 0 ) 

°-^-0.6 (1.5) 


vations will enable us to include a more exhaustive sample 
of Local Group galaxies. 


2.3 Uncertainties in, and limitations of, our 
approach 

There are a number of assumptions we must make in our 
modeling. For completeness, we list several of them here. 

• Star formation histories', we have only used the best-fit 
SFH, i.e., we assume perfect knowledge of the stellar mass 
at 2 = 5. This is mitigated by the SFHs reaching the ancient 
MSTO, but it still can introduce uncertainties at the factor 
of < 2 level into the total stellar mass formed. 

• Aperture corrections-, we assume the measured SFHs 
are representative of the entire galaxy. If the galaxies have 
strong population gradients, however, this assumption may 
be violated. The results of Hidalgo et al. (2013) indicate 
that at least some Local Group dwarf irregular galaxies do 
have population gradients, with younger stellar populations 
being more centrally concentrated than older populations. 
Any bias in our calculations is therefore likely to be in the 
direction of underestimating the UV luminosities at early 
times, as most of the galaxies have observations that sample 
close to their centers. 



Figure 3. Observed V-band magnitude at 2 = 0 versus modeled 
2 ~ 7 UV magnitude. Error bars indicate minimum 68 % confi¬ 
dence intervals (see Sec. 2.2 for more details). There is a well- 
defined, linear relation between the present-day V-band magni¬ 
tude and the modeled high-redshift UV magnitude for the galax¬ 
ies in our sample with a dispersion of ~ 0.5 magnitudes. The 
detailed form of this relation relies on the mode of star forma¬ 
tion in galaxies. For the dwarf galaxies studied in this paper, we 
have assumed that star formation occurs predominantly in short- 
duration bursts. More massive galaxies are likely to maintain star 
formation over longer periods, which would result in a shallower 
relation between Muv (2 ~ 7) and My (2 = 0). 


• Stellar IMF', we assume a fully populated Kroupa IMF, 
but the IMF may not be Kroupa or fully populated. Fur¬ 
thermore, we are basing the SFH on the IMF of lower-mass 
stars, but it was the high-mass stars in the early Universe 
that produced the UV light (Zaritsky et al. 2012). 

• Dust: we do not make any corrections for dust in our 
modeling. This is unlikely to be a poor assumption for early- 
time progenitors of the Local Group dwarf galaxies we are 
considering, as the vast majority are very metal-poor even 
today. More massive, and vigorously star-forming, galaxies 
at the epoch of reionization may be dusty (see, e.g., Finkel- 
stein et al. 2012; Watson et al. 2015). The galaxies studied 
here are expected to be low-metallicity and forming stars at 
rates of — 10“^M©yr“^ at 2 ~ 7, however, and are 
likely to be dust-poor (Fisher et al. 2014). 

• Metallicity: we assume a constant metallicity of 0.2^© 
throughout our analysis. Variations around this choice are 
unlikely to affect our modeling of Muv appreciably; see 
Johnson et al. 2013 for further details and a more expan¬ 
sive discussion of several of the uncertainties discussed here. 

• Distances: Our present analysis ignores the difference 
in luminosity distance over the width of photometric red- 
shift distribution P{z). For our 2 = 7 sample, the distance 
modulus only varies by ~ 0.2 magnitudes from 2 = 6.5 to 
2 = 7.5, meaning it will introduce at most a small correction 
to the effects we have modeled. The difference at 2 ~ 2 is 
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~ 0.4 magnitudes, which is still sub-dominant to the effects 
induced by bursty star formation. 

One potentially important uncertainty is the merger 
histories of dwarf galaxies. Our approach is inherently arche¬ 
ological in nature, as we are using resolved SFHs of galaxies 
observed at « = 0 to study the high-redshift Universe. Our 
analysis yields modeled UV fluxes as a function of time for 
the stars that are in the galaxy today, but those stars may 
have been formed in multiple distinct progenitors and only 
assembled more recently. If, for example, a typical merger 
history for our modeled galaxies were such that the 2 = 0 
galaxy had N equal-mass progenitors at 2 = 7, then our 
inferred UV luminosities at that epoch would be an overes¬ 
timate by a factor of N. We assess the likely distribution of 
progenitor masses in two ways, through the use of abundance 
matching to dark-matter-only simulations (for a statistical 
understanding) and through high-resolution hydrodynamic 
simulations (for individual case studies). 

Our first method relies on the ELVIS suite of dark- 
matter-only simulations of Local Group analogs (Garrison- 
Kimmel et al. 2014). We identify all objects within 1.2 Mpc 
of either Local Group giant at 2 = 0 (excluding the central 
subhalos, which correspond to the MW and M31) and trace 
these objects back in time. Rather than following only the 
most massive progenitor, we track all progenitors at each 
previous redshift. We find that, at 2 ~ 7, there is typically a 
dominant progenitor in terms of mass: the next most mas¬ 
sive progenitor is more than a factor of two lower in mass 
than the most massive progenitor at that epoch, on aver¬ 
age. Based on our modeling, this indicates that the main 
progenitor is likely to be at least one magnitude brighter in 
the UV than any other progenitor. Since the SFHs are highly 
episodic, this statement is true only in a time-averaged sense. 

We can also address this issue directly through hydrody- 
namical simulations. Specifically, we can use both an arche¬ 
ological and an instantaneous approach to 2 ~ 7 star forma¬ 
tion by asking (1) how many of the stars in the galaxy at 
2 = 0 formed by 2 = 7, and (2) what is the stellar content 
of the main progenitor of the 2 = 0 galaxy at 2 = 7? The 
difference between these two answers is a direct measure 
of the difference between resolved-star studies at 2 = 0 and 
observations at 2 ~ 7. We use the cosmological zoom-in sim¬ 
ulations described in Section 2.1 to perform this test. In all 
cases, the main progenitor of the simulated dwarf contains 
approximately 70% or more of the stars in the 2 = 0 galaxy 
that have formed by 2 = 7. The main progenitor is therefore 
already dominant by that time (see also Dominguez et al. 
2014), and our archeological approach to the properties of 
Local Group galaxy progenitors at high redshift should be 
a reasonable approximation. 


3 HIGH-REDSHIFT OBSERVATIONS IN A 
LOCAL GROUP CONTEXT 

With the results of Section 2, we are now in position to 
consider observations at high redshift in the context of well- 
observed galaxies in the Local Group. Of particular interest 
is the nature of the faintest objects observable in deep fields; 
in what follows, we discuss how the progenitors of Local 
Group dwarf galaxies would appear at 2 ~ 7 and 2 ~ 2. 


3.1 Local Group galaxies at cosmic dawn 

Figure 4 shows the observed UV LF of galaxies at 2 ~ 7 
(from Finkelstein et al. 2014; see also McLure et al. 2013; 
Schenker et al. 2013; Schmidt et al. 2014; Bouwens et al. 
2015b). The approximate halo masses expected for a given 
Muv are shown on the upper horizontal axis; these are ob¬ 
tained by matching cumulative number densities of galaxies 
and halos (assuming a Sheth et al. (2001) mass function 
for the latter). In the deepest blank field, the HUDF, HST 
is capable of reaching Muv ~ —17.5. Models of reioniza¬ 
tion in which galaxies play the dominant role In maintaining 
an ionized intergalactic medium require a major contribu¬ 
tion from fainter galaxies, however, meaning an important 
subset of the galaxies responsible for reionization have not 
yet been directly observed. The luminosities of the faintest 
galaxies contributing to reionization are unknown even the¬ 
oretically, as this limit (often denoted Mum) depends on the 
slope of the UV LF and the escape fraction of ionizing pho¬ 
tons from these galaxies (among other factors). The mini¬ 
mum Muv necessary for maintaining reionization in models 
ranges from Miim ~ —15 to Miim > —10 (e.g., Kuhlen & 
Faucher-Giguere 2012; Robertson et al. 2013). 

Our calculations of Muv ( 2 : ~ 7) for Local Group galax¬ 
ies can therefore provide crucial context for galaxies at cos¬ 
mic dawn. We find that the progenitor of the Large Mag¬ 
ellanic Cloud (LMC) - the brightest Local Group dwarf 
galaxy in our sample both at 2 = 0 and 2 ~ 7- had 
Muv(2: ~ 7) ~ —16, beyond the capabilities of HST in 
the HUDF^ The Small Magellanic Cloud (SMC) and IC 
1613 likely had Muv (2 = 7) ~ —13.8; virtually all reion¬ 
ization models therefore require progenitors of such galaxies 
to contribute to reionization. If Muv (2 ~ 7) « —13 galax¬ 
ies are required to maintain reionization, as is the case in 
many models, then the progenitors of the Sagittarius and 
Fornax dSphs must contribute. And finally, if galaxies as 
faint as M\jy{z ~ 7) ~ —9.5 are required, then progeni¬ 
tors galaxies such as Draco, CVn I, and Phoenix - which 
have Lv{z = 0) ~ (0.2 — 1) x 10® Lv,q - may be necessary 
contributors to reionization. 

These results also provide context for what observa¬ 
tions in the HST Frontier Fields and with JWST will see. 
At 2 ~ 7, Frontier Fields can probe galaxies as faint as 
Muv ~ —16.5, assuming magnifications of 5 in luminosity 
(1.75 mags); such observations may approach the depth re¬ 
quired for observing the main progenitor of the LMC. JWST 
deep field observations are expected to have a limiting mag¬ 
nitude of rriAB = 31.5 (Windhorst et al. 2006), correspond¬ 
ing to Muv ~ —15.5; JWST is therefore likely to reveal 
the 2 ~ 7 progenitors of Magellanic irregulars. If a Frontier 
Fields-like campaign with JWST could obtain a factor of 10 
in magnification (2.5 mags), it would observe objects as faint 
as Muv ~ —13; this would just approach the sensitivity re¬ 
quired to observe the progenitors of Fornax and Sagittarius, 


^ This prediction of Muv (2 ~ 7) ~ —16 for the LMC puts it 
very close to the UV luminosity for the MW, as predicted by 
B14 through abundance matching, at that time. This is consis¬ 
tent with expectations from the ELVIS suite (and other A-body 
simulations): the present-day most massive satellite is often only 
slightly lower (a factor of 1.5 — 2) in mass than the main halo 
itself at 2 = 7. 
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Figure 4. Local Group galaxies at 2 : ~ 7. The data points show the observed UV luminosity function at 2 ~ 7 from Finkelstein et al. 
(2014; HST, no lensing; black squares) and Atek et al. (2015; HST Frontier Fields; gray triangles). The best-fit Schechter (1976) function 
{M* = —21.03, 0* = 1.57 X 10“'* mag“* Mpc“®, a = —2.03) from Finkelstein et al. is displayed as a solid curve, while its extrapolation 
to lower luminosities is shown as a dashed line. The observational limits are also shown for the HUDF (assuming a completeness limit of 
Muv = —18), the HST Frontier Fields [assuming munj = 28.7 before lensing and a lensing magnification of 5 (1.75 magnitudes) or 10 (2.5 
magnitudes)], and JWST (assuming nin^ = 31.5; this same depth would be reached in the Frontier Fields at 2 = 7 with a magnification 
of 13.18 or 2.8 magnitudes). The modeled {Muy{z ~ 7)) for various Local Group galaxies are plotted as vertical bands (with error bars 
giving 1 O' uncertainties; see Table 1), indicating the power of local observations for interpreting deep-field data. Gonsensus reionization 
models require extrapolation to Muy{z ~ 7) ~ —13, corresponding to the brightest dSphs today (Fornax, Sagittarius). 


potentially revealing the faintest galaxies required for reion¬ 
ization. (We note that, based on the results of W14, we do 
not expect a strong truncation in the LF at at Muv ~ —13 
or even at significantly fainter magnitudes; the idea of a 
limiting magnitude Mum required for reionization is more of 
a mathematical construct than a physical cut-off. However, 
this does not preclude the possibility that the LF becomes 
shallower near Mum; see Sec. 3.1.) Progenitors of the vast 
majority of Local Group dwarfs will remain unobservable 
even with a JWST Frontier Fields-like project, however. 
This highlights the inherent difficulty of high -2 observations 
and the power of studying the high -2 Universe through its 
local descendants. 


3.2 Near-field / deep-field connections 

Figure 4 indicates that the faintest galaxies observable in 
the HUDF at 2 ~ 7 likely are hosted by Mvir « 3 x 10^° Mq 
halos, while the atomic cooling threshold of Tvir ~ 10* K 
corresponds to Muv ~ —10. B14 showed that the ELVIS 
suite of simulated Local Groups predicts approximately 50 
surviving, bound remnants of Mvir ( 2 : ~ 7) > 10® M© halos 
in the Milky Way’s virial volume today. They argued this 
was potentially problematic, as even low-level star formation 
in such halos would quickly over-produce the observed stellar 
content of Milky Way satellites. 

This tension is evident in Figure 5, which shows the 


2 ~ 7 UV LF of the Milky Way and its satellites (sym¬ 
bols) as well as predicted dark matter halo mass func¬ 
tions from the ELVIS simulation suite (gray shaded region). 
The corresponding values of Muv based on the abundance 
matching model described in Sec. 3.1 are given in the up¬ 
per horizontal axis. The LF from direct modeling of SFHs 
and from abundance matching are in good agreement for 
M^vy "L “12 (Mvir « 5 X 10® M©), but the disagreement 
disappears for fainter galaxies (lower mass halos), with low- 
mass halos far outnumbering the number of known galaxies 
even at the modeled 2 ~ 7 luminosity of Draco and Leo 
II (Muv ~ —9, corresponding to Mhaio ~ 3 x 10^ Mq). 
If every dark matter halo is capable of hosting a galaxy, 
then there should be 40-100 surviving descendants of galax¬ 
ies with Muv (2 ~ 7) < —10; our modeling predicts there 
are only 10 or so such galaxies around the Milky Way to¬ 
day. Either only a small fraction of the halos at this mass 
(Mhaio ~ 10® Mq) are capable of cooling gas and forming 
stars at 2 ~ 7 or the mapping between halo mass and UV 
luminosity is highly stochastic at early times in low-mass 
halos - both of which are contrary to current models and 
simulation results; or the UV LF breaks at Muv ~ —13, 
with Mvir < 10® Mq halos hosting fainter galaxies than our 
fiducial abundance matching model predicts. Whichever of 
these possibilities is correct, there are important implica¬ 
tions for the threshold of galaxy formation and the mass 
scale of halos that host classical and ultra-faint dSphs. 
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Figure 5. The z = 7 mass function of the main progenitors of 
surviving z = 0 (sub)halos - including the main progenitor of the 
MW itself - within 300 kpc of the Milky Way based on the ELVIS 
simulations (shaded region). The upper horizontal axis gives the 
abundance-matched A7uv; the data points show our modeled UV 
luminosity function for progenitors of the MW and its satellites 
at 2 ~ 7. The luminosity function from direct modeling of SFHs 
and from abundance matching are in good agreement for Mpv iS 

— 12.5 (Mvir ~ 7x 10® Mq) but diverge for fainter galaxies (lower 
mass halos), perhaps indicating the need for a break in the UV 
luminosity function at Muv ~ —13 at 2 ~ 7. 

The most realistic possibility may be that the 2 ~ 7 
UV LF breaks at Muv ~ ~13, as argued for in B14. Recent 
simulations have also found LFs that may flatten at fainter 
magnitudes than are probed by HST (Gnedin & Kaurov 
2014; O’Shea et al. 2015); there is no evidence for the LF 
flattening in current Frontier Fields data (Atek et al. 2015, 
Livermore et al. 2015, in preparation), which reach Muv ~ 

— 15.5 at 2 ~ 7. Muv ~ —13 is also a frequently-adopted 
value of Miim in models of reionization, although, as noted 
above, Mum need not be associated with any feature or break 
in the UV LF. 

If a break in the LF is indeed the answer, the Local 
Group data (in particular, galaxy counts at 2 = 0 com¬ 
bined with modeled UV luminosities at 2 ~ 7) indicate the 
faint-end slope for Muv ^ —13 should be close to — 1.2 
rather than the value of approximately —2 that is observed 
in the HUDF. As is shown in Figure 6, fainter galaxies would 
then live in halos that are more massive than our origi¬ 
nal abundance matching prescription would indicate. This 
issue is explored in more detail, and at a variety of red- 
shifts, in Grans et al. (in preparation). In this model (and 
in Fig. 6 ), the classical MW dSphs are hosted by halos with 
Mvir (2 ~ 7) = (0.5 — 1) X 10® Mq. All known Milky Way 
satellites could therefore be hosted by halos at or above the 
atomic cooling limit (Mvu ~ 10® Mq) at 2 ~ 7 (see also 
Mllosavljevic & Bromm 2014). 


Figure 6. Similar to Figure 5, but assumes a UV luminosity 
function that breaks to o = —1.2 at Muv > —13 (from the 
fiducial value of o = —2.03 for brighter galaxies). The 2 = 7 
census of galaxies surviving to 2 = 0 in the Milky Way is in 
much better agreement with the modeled UV luminosities in this 
case, as galaxies over a wider range in luminosity are placed in 
a narrower range of halo masses (for Muv > —13). In such a 
scenario, all known MW galaxies, including ultra-faint dwarfs, 
could lie at or above the atomic cooling threshold of Mvir ^ 
10® Mq at 2 = 7. 

It is important to emphasize that completeness in the 
2 = 0 data is not an issue when constructing Figures 5-6: we 
have only used data for satellites with Ly{z = 0) > 10® Lv,o 
and current Galactocentric distances of < 300 kpc, a re¬ 
gion where our census of satellites is very likely complete 
(e.g., Tollerud et al. 2008; Koposov et al. 2008; Walsh et al. 
2009). Furthermore, the mismatch in Figure 5 is already 
significant for galaxies such as Leo I and Sculptor (with 
Lv{z = 0) « 5 X 10®Lv,q). The only galaxy this bright 
that has been found within the Milky Way’s virial volume 
since the 1950s is Sagittarius (Ibata et al. 1995), whose pres¬ 
ence had been concealed by the Galaxy’s disk. Although we 
do not have data for Sextans, the discrepancy in numbers 
shown in Figure 5 is an order of magnitude, indicating that 
even the inclusion of Sextans and the discovery of several 
10 ® Lv ,0 satellites would not change the qualitative picture 
described here. 

3.3 Local Group galaxies at cosmic noon 

Dwarf galaxies are also important test-beds of galaxy for¬ 
mation physics at eras other than the epoch of reionization. 
In particular, large samples of UV-selected galaxies at 2 ~ 2 
are becoming available, and ongoing HST programs are en¬ 
abling studies of galaxies that are intrinsically as faint as 
Muv ~ —13.5 (Alavi et al. 2014) at that time, which is 
close to “cosmic noon”, the peak of the cosmic star forma- 
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Figure 7. Local Group galaxies in a z ~ 2 context. The data points show HST observations of z ~ 2 galaxies in the GOODS-South field 
(Oesch et al. 2010; black points) and in the A1689 lensing field (Alavi et al. 2014; gray points). The best-fitting Schechter luminosity 
function from Oesch et al. (M* = —20.16, <f>* = 2.188 X 10“^ mag“^ Mpc“^, a = —1.60) is plotted as a solid curve, and its extrapolation 
to fainter magnitudes is plotted a dashed curve. Vertical bands indicate the modeled {Muy{z ^ 2)) for various Local Group galaxies (with 
error bars giving 1 cr uncertainties; see Table 1). While deep observations of blank fields only probe much more massive galaxies than 
those of the Local Group (except the Galaxy and M31), lensing magnification likely enables the study of galaxies similar to progenitors 
of the LMG and SMG at z ~ 2 and are nearly capable of reaching the main progenitor of galaxies like IC 1613 and Fornax at that time. 


tion rate density (Madau & Dickinson 2014). Understanding 
the likely descendants of such galaxies today - or the likely 
progenitors of Local Gronp galaxies - will shed farther light 
on the processes at work in galaxy formation over the past 
10 billion years. 

Fignre 7 shows the observed UV LF of galaxies at 
z ~ 2 (black data points, from Oesch et al. 2010). The 
gray data points are observations that take advantage of 
lensing magnification combined with deep near-UV imag¬ 
ing (WFC3/UVIS in F275W), which allowed Alavi et al. to 
probe to much fainter galaxies (Muv ~ ~13) than would 
otherwise be possible. The figure also shows the UV magni¬ 
tudes that a variety of Local Group galaxies would have at 
this redshift. 

The LMG and SMC are predicted to have M\jy{z ~ 
2) « —16 and —14, respectively, placing them well within 
the reach of HST at z ~ 2. IC 1613 and the Fornax dSph 
are predicted to have Mijy(z ~ 2) « 12.5, meaning the 
progenitors of these galaxies are just at the edge of HST’s 
current capabilities at z = 2. Other Local Group galaxies in 
our sample are all predicted to be fainter than Muv = —12 
at z ~ 2 . 

As yet larger samples of galaxies become available, it 
will be possible to examine z ~ 2 galaxies as a popula¬ 
tion in ever-greater detail. If lensing allows the study of 
Muv ~ “12 galaxies at z ~ 2 , a number of applications 
will present themselves. Such observations will likely reveal 
the progenitors of present-day dwarf irregular galaxies and 
massive dSphs in a variety of environments. The distribution 
of z ~ 2 UV fluxes obtained through direct observation will 
help constrain potential burst parameters (e.g., the z ~ 2 


equivalent of Figure 2). Observations at z ~ 2 will also in¬ 
crease the precision with which the faint end slope of the 
UV luminosity function can be measured, placing further 
constraints on the evolution of cosmic star formation and 
on the properties of dust in faint, star-forming galaxies. 


4 DISCUSSION 

The previous sections illustrate the power that deep 
resolved-star observations of Local Group galaxies have for 
understanding not just the nearby Universe but also (much) 
earlier eras. Using the modeled UV luminosities, we have 
shown that Local Group galaxies such as the Fornax dSph 
and IC 1613 were, at z ~ 2 , slightly less luminous (Muv ~ 
— 12.5) than the faintest galaxies HST is currently observ¬ 
ing at that epoch (Muv ~ —13). The LMC and SMC are 
predicted to be substantially more luminous at that time; 
many of the faint (Muv ~ “15) galaxies observed at z ~ 2 
may evolve to become Magellanic irregulars in the local Uni¬ 
verse. Pushing to higher redshifts (z ~ 7), we find that 
none of the Local Group galaxies are bright enough to be 
seen in the HUDF, consistent with expectations based on 
number counts and space densities (e.g., Trenti et al. 2010, 
B14). Many reionization models predict that extrapolating 
the observed LF to Muv ~ “10 is required to maintain 
reionization; if this is the case, the progenitors of some of 
the faintest classical dSphs (e.g., Draco, CVn I) and low- 
luminosity dwarf irregulars (e.g., Phoenix and Leo A) would 
be required contributors to the ionizing background. 

In Section 3.2, we showed that a comparison between 
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our derived Muv values for Local Group galaxies at z ~ 7 
and counts of surviving subhalos in dark-matter-only sim¬ 
ulations of the Local Group show that 2 ~ 7 LF cannot 
continue with a slope of a = —2 to Muv ~ —9, as this 
would imply an order-of-magnitude excess in surviving satel¬ 
lites with z = 0 luminosities comparable to the classical 
MW dwarfs. The faint-end slope derived from high-z obser¬ 
vations is still somewhat degenerate with measured value 
of M* in the Schechter LF, however, with a typical uncer¬ 
tainty of ±0.3 (e.g., McLure et al. 2013; Finkelstein et al. 
2014; Bouwens et al. 2015b). These degeneracies translate 
into substantial uncertainties in Mum (see, e.g., figure 2 of 
Robertson et al. 2013), which in turn will affect the specific 
value of Muv at which archeological studies in the MW re¬ 
quire the high -2 LF to break. The general trend is that if 
the faint-end slope of the LF is steeper, it must break at 
brighter values of Muv; if it is shallower, the break can oc¬ 
cur at fainter values. This mirrors the behavior of Mum and 
its dependence on a. 

The Planck collaboration has recently reported a value 
of the optical depth to electron scattering, r = 0.066 ±0.012 
(Planck Gollaboration 2015), that is « 30% lower than 
previous determinations of r = 0.088 ± 0.014 (Komatsu 
et al. 2011). Robertson et al. (2015) show that extrapolat¬ 
ing the LF to a limiting magnitude of Miim « —13 (with 
a luminosity-independent escape fraction /esc = 0.2) is suf¬ 
ficient to maintain reionization (see also Finkelstein et al. 
2012) and match the updated Planck determination of r. 
Such a scenario still requires extrapolation of the LF down 
to progenitors of the Fornax dSph, which has a present-day 
luminosity of Lv ~ 10 ^ Lv,q- 

If time-averaged escape fractions are lower than 20% for 
galaxies with Muv ^ ~13, as suggested by the work of Wise 
et al. (2014) and Ma et al. (2015), then maintaining reion¬ 
ization would require even fainter galaxies (or additional 
sources such as X-rays; e.g., Mirabel et al. 2011). Galaxy- 
driven reionization scenarios therefore still require that most 
of the known Local Group irregular galaxies - and the most 
luminous of the dSphs - were necessary contributors to cos¬ 
mic reionization. These galaxies are intrinsically very faint 
at high redshift: JWST will only resolve galaxies that are an 
order-of-magnitude brighter even in a deep field campaign. 
The unique high-angular-resolution capabilities of observa¬ 
tories such as HST and JWST, and the deep observations 
of faint galaxies in and around the Local Group they facili¬ 
tate, therefore hold the promise of providing a unique probe 
of the earliest epoch of galaxy formation. 
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